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Abstract 

We present a new derivation of the Berry phase picked up during ex- 
change of parallel vortices. This derivation is based on the Bogolubov - 
de Gennes formalism. The origin of the Magnus force is also critically 
reanalised. The Magnus force can be interpreted as interaction with effec- 
tive magnetic field. The efective magnetic field may be even of the order 
10 6 T/A. We discuss a possibility of the FQHE in vortex systems. As the 
real magnetic field is varied to drive changes in vortex density, the vortex 
density will prefer to stay at some quantised values. The mere existence 
of the FQHE does not depend on vortex quantumstatistics although the 
pattern of the plateaux does. We also discuss how the density of anyonic 
vortices can lower the effective strengh of the Magnus force, what might 
be observable in measurements of Hall resitivity. 

To appear in Physical Review B. 
cond-mat/9508140 



1 Introduction 

In our recent paper 0] we have pointed out the possibility that vortices in super- 
conducting films might be anyons. There is a classic paper by Haldane and Wu 
[0 , which demonstrates that vortices in superfluid helium layers are not anyons 
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because there is no well denned dilute limit. We have reanalysed this ques- 
tion from the point of view of the phenomenological time-dependent Ginzburg- 
Landau models. It was shown that there are situations when one can go around 
the argument in Q. A classic example are vortices in the Ginzburg-Landau 
models of the fractional quantum Hall effect j^, |j. Chern-Simons interaction 
makes them well localised objects and the dilute limit can be defined. The main 
effect of the CS gauge field is to remove the divergent gradient energies present 
in the global model. These divegencies are also removed in the Ginzburg-Landau 
model for superconductors. We have considered a gauged nonlinear Schrodinger 
equation as a minimal version of time-dependent Ginzburg-Landau model. To 
model the structure of the vortex core we assumed it was filled with normal 
fluid so as to make the whole structure locally charge-neutral. The strenght of 
the statistical interaction was proportional to the net deficit of the superfluid 
replaced inside the core by the normal phase. 

It is argued || that the time evolution of the condensate is much faster then 
that of the normal fluid. Thus in the static case the core is filled with the normal 
fluid but when the vortex moves fast enough the normal fluid can not adjust 
itself and does not follow the vortex motion. It is a crude approximation |^| and 
it does not invalidate our results. The statistical interaction shows up in the 
adiabatic approximation which is quite an oposite limit. In this limit we can 
assume the nonuniformity of the normal fluid does follow the vortex motion. The 
effective Lagrangian for vortex motion can be arranged term by term according 
to the powers of vortex velocity. The statistical interaction together with the 
term responsible for the Magnus force Q are the first terms in this effective 
Lagrangian. The distortion of the normal fluid distribution from the static one, 
which is at least of the first order in velocities, can contribute but only to the 
term quadratic in velocities which is the next to leading term. Thus at least for 
slow motion as compared to characteristic velocities, when expansion in powers 
of velocities is justified, the Magnus force and statistical interaction are not 
altered. 

One can rewrite the local BCS model in terms of the gap-function field 
[H 0- P but the effective theory appears to be nonlocal - it contains derivatives 
of arbitrary order. In this way one goes from description in terms of electronic 
degrees of freedom to description in terms of Cooper pairs. The latter is well 
suited for the bulk of the superconductor. However in the core of the vortex we 
can expect some decoherence, which phenomenologically might be decribed by 
normal fluid. Whenever degrees of freedom other than those of Cooper pairs 
come into play the description in terms of Cooper pairs may happen to be 
irrelevant. One way of dealing with the problem is to truncate the effective 
theory on the lowest order in derivatives and introduce more or less explicitly 
something like a normal component. This approach is limited by the poor 
knowledge about the nature of the normal fluid. Another approach is to take 
as many orders in derivatives in the effective theory as possible to describe also 
the normal fluid in terms of the gap function. The problem is that in practice 



2 



one would have to cut the expansion at certain order and there is no warranty 
that such a cut theory would be self-consistent. 

To go around these problems we will derive the statistical interaction in 
the framework of the Bogolubov-de Gennes formalism for pure samples at zero 
temperature. We will work in the quasi two-dimensional regime of long parallel 
vortices. In the case of superconducting layers thicker than 100A the penetration 
lenght A is still very close to the penetration lenght A in the bulk superconductor. 
In this regime we expect modifications to be rather quantitative in nature then 
qualitative. An important first step was done by Gaitan [fl0[ in his derivation 
of the Berry phase responsible for the Magnus force. In this paper we are going 
to reanalize his derivation and then to extend the method to the case of two 
well separated vortices. In the microscopic theory the language of the normal 
and superfluid is not very fruitful. The distinction of the states below the Fermi 
surface into localised bound states and scattering states appears to be more 
natural. The distinction does not influence the value of the Magnus force but 
it is crucial for the statistical interaction. The scattering states are common to 
all vortices while bound states can be identified with particular ones. 



2 Preliminaries on vortex solution in BCS the- 
ory 

The problem of the vortex solution in the BCS theory can be conveniently posed 
within the Bogolubov-de Gennes formalism p"2| . It has not been completely 
solved although some qualitative fictures of the solution are known pl| . We 
will restrict here to listing the basic ingredients of the formalism. 

The Bogolubov equation, defined in the Nambu spinor space is 
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where Hbog is the Bogolubov hamiltonian 



Hbog = 



-i(V - ieAf - E F 



A(x) 
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Ep is the Fermi energy. A(cc) = Ao(r) exp(— iO) denotes the gap function of 
vortex solution. A (r) interpolates between at the origin and a constant, 
which we will call -^po, at infinity. There are both positive and negative energy 
solutions. If (u n ,v n ) is an eigenstate with energy E n > 0, then (— v*,u^) is a 
solution with energy —E n < 0. The equations (|lj) have to be supplemented by 
a self-consistency condition 



(3) 
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where g is the BCS coupling constant, together with Maxwell equations deter- 
mining the vector potential. 

The field operator for Nambu quasiparticles can be expanded in terms of the 
solutions of Eq. ([!]) . 



•<•>-( $5 )-£ 
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(4) 



The negative energy states are occupied in the BCS ground state 

\BCS>=Y[ 7nl \0> . (5) 

n 

The eigenstates satisfy the following orthogonality 

tfx [u n {x)u* m {x) + v n {x)v* m (x)] = S nm , 

d 3 x [u n {x)v m {x) - v„(x)u m (x)} = (6) 

and completeness relations 

^2[u n (x)u*(x') + v^(x)v n (x')] = S(x - x') , 

71 

J2[u n (x)v*(x') - v*(x)v n (x')} = . (7) 

n 

With the help of these relations the creation and annihilation operators can be 
expressed as 



Inl = J d 3 X [-^\(x)vl{x) + i>l(x)u*(x)] , 

lU = / d3x &\(n)<(v) + i>l(x)Vn{v)] ■ (8) 



Adiabatic vortex motion jlO| gives rise to a Berry phase in the solutions of 
Eq-(jlj), (u n ,v n ) — + exp[i<f> n ](u n , v n ). These Berry phases sum up to the total 
Berry phase picked up by the ground state | BCS >^ exp[ir] | BCS > which 
with the help of Eqs.(|||5|) can be established to be 

r = -£>„ . (9) 

n 

To persue some of the questions we need more detailed knowledge about the 
eigenstates. The axially symmetric ansatz takes the form 

Xn(x) = ( U v j^ ) = e^-e^-^Ur) , (10) 
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where Tik z is the z-component of the momentum, fi must be a half-integer for 
the expression to be single- valued. The functions f n (r) have been investigated 
in with the help of WKB approximation. We will quote some more detailed 
results in the following. 



3 Origin of the Magnus force 

Now we are going to rederive the Berry phase responsible for the Magnus force 
following the argument of Gaitan [[To) . In comparison with jl0| we clarify some 
points and remove some unnecessary assumptions. 

The general form of the Berry phase in two dimensions is 

(t> n = ijdt J d 2 xxl(-^+t^A )xn ■ (11) 

The time derivative is understood as a total derivative with respect to slow 
degrees of freedom. For an adiabatic motion of a single vortex the derivative 
has to be replaced by foVr , where t*o is a position of the vortex singularity. 
The scalar potential vanishes for the vortex solution so we will skip the second 
term in what follows. With the axially symmetric ansatz (|l0|) the phase becomes 

4> n = Jdt J d 2 x r„[0d(-/i + \<T z )Xn)V ro + fl(r)Vr f n (r)} (12) 

The second term vanishes by symmetry arguments. The contribution of the 
first term to the total Berry phase is 

r = - E = - I dt [ d 2 x (r o V ro 0)S(r) (13) 

n J J 

where 

%)=E0^W! 2 (-M+5)+l^(r-)| 2 (-M-5)] • (14) 

n 

TiS is minus the z-component of the canonical angular momentum density. It is 
not a gauge-invariant integral of motion. S is the expectation value density of 
the operator —ih-^g instead of the gauge-invariant —ihJjg + o z Aq. 

S(0) = because either u n (v n ) or the factor (— /i — |) vanishes at the 
origin. To find out its asymptotic behavior at infinity we would need a much 
more detailed knowledge about the solutions. We go around this problem by 
resorting to the effective theory which is equivalent to the microscopic formalism. 
The general term linear in the covariant time derivative reads 

Jdt Jd 2 Xim W t + ^A )A-c,.}[G(A*A) +sp aUalder iV aU V eterm S ] . 

(15) 
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By "spatial derivative terms" we mean terms which are of at least first order 
in the covariant spatial derivatives. G is a function of A* A only which tends 
to Ps/po as the gap function A approaches its asymptotic equilibrium value. p s 
is the equilibrium Cooper pairs' density. Let us consider the adiabatic rotation 
of the vortex solution A = Ao(r)e"' 9 around its axis, 9^9 — tot. The action 
picks up a term (to lowest order in to) 



,JdtJ { -RA*A[G(A*A) + spatial derive terms]} 



(16) 



The spatial integral is just the total angular momentum. For large r, where Aq 
tends to zero, the density of this angular momentum is, by gauge invariancc 
(151), equal to h times minus the bulk Cooper pairs' density HS « — hp s = 



h]hn r - >00 ^eAo) ' w h ere W 1S the effective action and p s > 0. We have to 
stress that we make use of the effective theory only very far from the vortex 
core where it should be equivalent to the microscopic treatement. In particular 
in the distant asymptotic region there is no contribution from unpaired bound 
states which can not be described in terms of Cooper pairs. 

We have all we need to calculate the Berry phase. Let us expand the in- 
tegrand in Eq.(|l3|) around the vortex position ro = (X,Y) close to the origin, 

(x,y) = (o,o), 

S = S(r) - S'(r)[X cos9 + Y sin 6} + 0(rg) , 

■ „ „ v,>n6> X sin 29 - Y cos 29 . ■, cosO X cos 29 + Y sin 29 . 

r V ro 9 = X(— + ) +Y{-— ^ ) + 0(rg)7) 

A straightforward integration yields 

F = -7r[5(oo) - S(0)} f dt e kl X k X l + O(r 2 ) . (18) 



The expression O(rp) does vanish. We are considering single vortex in absence 
of any driven current. Such a system is translationally invariant and isotropic. 
The first term on the R.H.S. of Eq.(|lJ) is already the most general term linear 
in velocity which is, up to a total time derivative, translationally invariant and 
isotropic. Thus we do not need to consider finite ro to obtain a generally valid 
expression. The phase (|l|) is remarkably simple to evaluate. For a vortex with 
winding number —1 it reads 



T = tt Ps ] dte kl X k X l . (19) 

From our derivation of the Magnus force it is clear that the Wess-Zumino term 
(in the gauge A — 0) 

J dt J d 2 x[ Ps d t 9] , (20) 



G 



with p s = const, does not make much sense as it stands. p s is the same at 
the origin as at infinity so the Magnus force vanishes. The formula (j2TJ) is to 
be understood with an implicit assumption that a small area around the phase 
singularity is excluded from the spatial integration. In other words p s must be 
put equal to in this area. It is not difficult to realise, by performing radial 
integration first and then integration over the angle around the singularity, that 
the way of regularisation does not matter. In particular it does not need to be 
rotationally symmetric. The only factors that determine the Magnus force are 
the two limit values of p s . Thus vortices in a condensate will always feel the 
Magnus force. It is not the case for say Jackiw-Pi solitons fl(| , where p s is zero 
both at the origin and at infinity. 



4 Mutual statistical interaction of vortices 

Let us consider two vortices: "1" at the origin and "2" very far apart at R(t). 
It is important to realise that the eigenstates of the Bogolubov hamiltonian 
can be divided into common scattering states, which we will still denote by just 
u n , v n , and bound states which can be identified with a given vortex Un , Vn • 
Vortices are very distant so there is no overlap between their localised bound 
states. 

The bound states of the stationary vortex " 1" feel what is going on around 
them through the pair potential A(t,x) inside and around the core. Vortices 
are well localised so a fairly good approximation to a two-vortex gap function 
is the product ansatz 

^A(t,x) = A v (x)A v [x-R(t)} , (21) 

where A v (x) = A (r) exp(— iO) denotes the gap function of a single vortex 
centered at the origin. Close to x = this expression can be further simplified 

A(t, x) = A v {x)e~ l9[x ~ R(t)] , (22) 
Thus the bound states of the static vortex have to be modified as 

X% ] [x,R{t)] = e~^ x ~ R(t)] X n{x) . (23) 
Their contribution to the Berry phase is 

-Jdt Jd 2 x {RV R 9[x - R(t)]} ]T {\ | u n | 2 - 1 - I v n | 2 ) w 

bound st. 

S \\ Sx Ed M « I" - I w « I')} / dt &V R 9[R(t)] , (24) 

bound st. 

where the approximate equality is valid for small % , where r c is a radius of the 
core. The equality | u n \ 2 — \ v n | 2 holds for the bound states, at least up to the 
WKB approximation urn, so their contribution to the Berry phase vanishes. 
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Now as the vortex "2" moves its bound states follow its trajectory x^ [r — R(t)], 
similarly as in the single vortex case considered in the previous section. In ad- 
dition, as an effect due to the vortex "1" (pi]), their components perform the 
relative phase rotation 

X ^[x,R(t)]^e-^ R ^ Xn [x-R(t)} . (25) 

The contribution from this relative phase rotation is once again zero. The bound 
states contribute but only to the Magnus term (|lj]) just as in the single vortex 
case. 

The bound states do not give rise to any new effects so let us consider scat- 
tering states common to both vortices. In the vortex core region the asymptotes 
of the scattering states must be close to those of the scattering states for a single 
vortex but the phase has to be replaced by the asymptote of the phase in the 
product ansatz (p2|). Close to the origin 

Xn [x, R(t)} « e *.* e i<J>- s *)MB)+0l* ~ R Wf n (r) . (26) 
The contribution from around the stationary vortex is 

- / dt J d 2 x KV R e[x - R{t)]S{r) , (27) 

where S is a part of the canonical angular momentum due to the scattering 
states 

S(r)= l\Mr)\ 2 (-V+\)+\v n (r)\ 2 (-»-\)] . (28) 

scatt. st. 

Far from the core S ~ S ~ —p s - If S were equal to — p s also in the core, the 
contribution to the Berry phase from (|27j) would be just the same as to the 
Magnus term. Thus we are interested only in the effects due to deviations of S 
from its asymptotic value — p s . Inside the core | u n | 2 and | v n | 2 are changed by 
a factor which is > 1 for the states with [i negative and < 1 for /i positive . 
The net deviation SS(r) = S(r) + p s is positive. At the very origin 5S(0) = p s . 
The total change in the Berry phase is twice that in (|2^), as there are two 
vortices, and amounts to 

ST = [2 J d 2 x SS(r)} J dt RV R 6{R) . (29) 

Thus the total Berry phase for a dilute vortex system is 

f d 

r = / dt l-n Ps n p e a X k (p) X\ p) + a ]T n p n q -G {p , q) ] , (30) 

p p<g 

where the indices p, q run over vortices, n's are their winding numbers, 6( p , 9 ) is 
the angle between the p-th and q-th vortex and the numerical factor a can be 
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read from Eq.(]29|) 



a = 2 / d 2 x 5S{r) > . 



(31) 



a is roughly the number of electrons inside the core and as such it can range 
from ~ 1 for high T c superconductors to ~ 10 5 for some conventional type II 
superconductors . 

5 Vortex statistics within variational wave- function 
approach 

Once we have derived statistical interaction in the microscopic setting it may 
be worthwhile to reanalise some earlier approaches to similar problems. In the 
paper by Ao and Thouless j?]] the Magnus force was derived with the help of 
the variational many-electron vortex wave-function 



where zq is a complex vortex position, z^s are positions of electrons and ^po[z\ is 
an antisymmetric variational function. The phase factors in the wave-function 
are determined by the demand of correct electronic quantumstatistics and by 
topological properties. There are variational profile functions in ip which can 
not be established without dynamical considerations. One can consider an adia- 
batic vortex motion along some trajectory and calculate the Berry phase picked 
up by the wave-function. This Berry phase coinsides with Eq. ( fl9| ) . For a closed 
path the Berry phase is proportional to the number of electrons enclosed by the 
trajectory. 

This setting is convenient to analise what is the dependence of the Magnus 
force on impurities 0. An impurity can be viewed as an attractive potential 
which traps some of electrons in localised bound states. The trapped electrons 
disappear from the ansatz ( |32"|) . The Berry phase is still proportional to the area 
enclosed by the trajectory but this time the area should not be multiplied by 
the total density of electrons but rather by the total density minus the density 
of electrons trapped by impurities. Impurities lower the value of the Magnus 
force. 

Now let us consider the effect of an exchange of two vortices. More precisely, 
let us fix the position of one vortex and consider another distant vortex moving 
around it. One could argue there is no special effect because the net charge of 
any vortex must be zero. Provided the trajectory is large enough, there is no 
change in the number of enclosed electrons due to the enclosed vortex. The last 
sentence is certainly true but the example with impurities tought us that it is 
not the total number of electrons that really matters but rather the number of 
electrons in the coherent state described by the wave- function (p2|). We know 




(32) 



k 
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from the discussion in the previous sections that inside vortex core the scattering 
or continuum states are replaced by bound states. Thus vortex can be viewed 
as a kind of impurity, which traps some of electrons into localised bound states 
with energies within the energy gap band. The localised electrons are removed 
from the wave-function (^). There is an additional Berry phase proportional 
to the number of vortices enclosed by the trajectory. Each enclosed vortex 
contributes a term proportional to the number of electrons trapped inside its 
core. 

We can consider a path for a chosen vortex in a more or less uniform dis- 
tribution of vortices. If we neglect possible intervortex correlation effects, the 
background vortices could be regarded as uniform distribution of impurities 
lowering the density of electrons in the coherent state (|32|). In this mean-field 
approximation the Magnus force acting on a choosen vortex is lowered by the 
presence of another vortices. This approximation is nothing else but the delo- 
calisation procedure so often applied to anyonic systems. The Magnus force can 
be interpreted as Lorenz force due to interaction of effectively charged vortices 
with some uniform effective magnetic field. The statistical interaction can be 
seen as Aharonov-Bohm effect due to the fluxes attached to vortices. In the 
mean-field approximation the fluxes, which are opposite to the external flux, 
are delocalised and they lower the net uniform flux. In the same way the real 
impurities can be interpreted as localised fluxes, opposite to the external field, 
randomly distributed over the plane. If the M-F approximation appears to work 
for real impurities, it will also work for vortices. 



6 Hall angle and vortex density 

The fact that the value of the Magnus force can be lowered with increasing 
density of vortices can, in principle, be observable in Hall experiments [ pT[ . The 
vortex density should increase and the M-F Magnus force should decrease with 
increasing real magnetic field. This should manifest itself in the changes of the 



measured Hall angle. Vortex equation of motion takes the form 15 

m e //f = ^2~(r — v a ) X z — rjdr + F pin + f , (33) 

where m e // is a small effective vortex mass, 77 is a vortex viscosity, F p i n is 
a pinning force, / is a fluctuating force, d is a sample thickness and v s is a 
driven uniform superfluid velocity. When we neglect pinning and average over 
fluctuations the stationary state motion will be determined by the equation 

(f - v a ) X z = t'An(6 H )r , (34) 

where tan(##) = -=^. The solution is 

iU _ v s + (zx v s )tan{0 H ) , » 

r l + tan 2 (M ' [ ' 
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Oh is the angle between the superfluid velocity v s and the stationary vortex 
velocity. The angle is the larger the weaker is the Magnus force. If the effective 
Magnus force is lowered with increased vortex density the angle should also 
grow with external magnetic field, which drives the rise in vortex density. The 
changes of the Magnus force due to changes in vortex density should be the 
more rapid the larger is the number of electrons trapped in the vortex core. For 
this reason we would recommend experiments on mildly type II conventional 
superconductors with large vortex cores (large correlation length £). Rather 
strong viscosity should be prefered for the angle to be more sensitive to the 
strengh of the Magnus force. The sample should be pure of pinning centers to 
avoid obscure pinning effects. 



7 Vortices' fractional quantum Hall effect 

To summarise our knowledge about the dynamics of planar vortices, let us write 
down an effective Lagrangian for diluted vortices with topological charge —1 

p p<q 

(36) 

The indices p, q run over vortices. m e f f is the effective vortex mass. It is usually 
estimated to be around 10 8 m e /m = 10~ 2 m e / A. 

The second term in Eq.(^6|) decribes interaction of vortices with effective 
magnetic field. We stress that this field has nothing to do with the real magnetic 
field B ext , which in this case is just a device to drive the changes of vortex 
density. If we assumed the density of electrons to be ~ 10 30 rn~ 3 = lit -3 the 
effective magnetic field defined by eB ^ ,f = irhp s would turn out to be ~ 10 6 T/A. 
When compared with the effective vortex mass per lA the magnetic field turns 
out to be incredibly strong. Its effect on a vortex should be the same as that of 
the 10 8 T magnetic field on an electron. Vortices can be expected to be confined 
to the lowest Landau level (LLL). 

Now let us consider a single vortex at z = X\ + 1X2- What is the magnetic 
lenght I which determines the size of the LLL wavepacket 

^"-vB?*"*' (37) 

The magnetic lenght is determined by the strength of the exactly known Magnus 
force I 2 = 7t^— ■ The area over which the center of such a vortex fluctuates can 

2-irps 

be estimated to be — , which is of the same order as the area per one electron. 
This effect is significant for extremely type 77 superconductors where the core 
is very thin. 

Now we are prepared to address the question of the fractional Hall effect. 
Vortices are anyons with a statistical parameter a. A nontrivial statistics would 
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be sufficient to prevent them from overlapping if the third term in ( pq ) were not 
regularised at short distances and replaced by mutual charge-Mux interaction 
[ p~7[ . Fortunately we also have an effective short range mutual repulsion V e fj 
which at low temperatures may be sufficient to keep vortices at a distance. On 
the other hand the potential is very weak as compared to the Landau energy so 
the mixing with higher LL's should be in any case negligible. Following Laughlin 
[ fl8| , the trial wave- function for a many vortex state can be written 

4> m [z] = l[(z p - z q ) a+2m i[exp(- | z r | 2 /Al 2 ) , (38) 

p<q r 

where positions of vortices, m is a nonnegative integer so that the 

exponent (a + 2m) provides a correct quantumstatistics. The density n and the 
filling factor v in such a state are 

v m = 2irl 2 n m = — . (39) 

For m — the density is just no — 1/2-Kal 2 — p s /a. a is roughly equal to 
the number of electrons inside the core so in the m = state vortex cores 
would have to overlap slightly. Certainly this density is close to that of the 
Abrikosov lattice E^] . For larger m the density is smaller and finally we should 
get outside of the crystalline regime. Then as the density of vortices is driven to 
change by the changes of the real magnetic field B ext we should observe some 
plateaux at the densities n m and maybe also at some other quantised filling 
factors. For conventional superconductors, which are mildly type II, a is large, 
so the difference between n m and n rn+ i is too small for the plateaux to be 
distinguishable experimentally. For extremaly type II superconductors or high 
T c superconductors we can expect a to be even ~ 1. In such a case the first 
plateaux can be expected already at ~ 1/3 and ~ 1/5 of the Abrikosov lattice 
density. 

The main observation of this section is that if the Magnus force is translated 
into the language of interaction with some effective uniform magnetic field, the 
field appears to be even ~ 10 6 T/A. It has nothing to do with the real magnetic 
field which is at best ~ 10T. The FQHE is a result of this huge effective magnetic 
field and of repulsive intervortex potential but its existence does not depend on 
the vortex statistics as the FQHE is possible even for bosons [^0| . However the 
pattern of the FQHE plateaux can help to identify the quantumstatistics. 

An experimental setup to detect the FQHE would consist of a planar sam- 
ple of some pure superconducting material in external uniform perpendicular 
magnetic field. One would have to measure the total flux which penetrates 
through the sample. This flux gives the actual number of vortices in the sample 
because each vortex carries one flux quantum. Just after continuously turn- 
ing on weak magnetic field the field lines would be pushed out of the sample. 
Only above some threshold value of the flux the energy of the field could be 
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minimised by creating a vortex line. The story would repeat until some stable 
FQHE plateau were reached. The plateau is a manifestation of a very stable 
many-vortex state so an addition of one vortex to this state would cost more 
energy then to a state far from the plateau. When the stable state is approached 
from low densities the addition of one more vortex should be much easier than 
usually. Quite opposite, when we lower the external magnetic field it is easier 
to remove one flux quantum from the sample just above the stable state and 
more difficult to remove it just below the state, see the figure. The moduli of 
the change in the driving flux A<£> (or some equivalent quantity) necessary to 
change the flux through the sample by one quantum will develop a characteristic 
pattern around the stable density no. Two measurements, one with adiabati- 
cally increasing and one with adiabatically decreasing external magnetic field, 
would give two curves with opposite polarisation. Taking their difference will 
amplify the effects due to the plateau and remove the not neccesarily constant 
bias. 

It should be stressed here that our results are not in contradiction with pre- 
dictions of a bosonic Hall effect in Josephson junction arrays ^l). The arrays 
are strictly planar devices. The penetration lenght A is likely to be greater then 
the sample size. What is more, it seems to be possible to excite a vortex with- 
out any bound states in the core, which indeed might be a boson. Our results 
are exact in the limit of long straight-linear parallel vortices. In practice it is 
sufficient that the sample is thick enough for the penetration lenght A to be 
close to that in the bulk. The sample thickness would have to be a bit more 
then d — 100 A. For such a d there is still no space for vortex entanglement. In 
the absence of entanglement vortices can be uniquely projected on a plane. In 
nonzero temperature one should expect transversal modes to be excited. These 
excitations do not affect the Berry phase picked up during vortices exchange as 
it depends on the overall topological properties of vortices. Vortices excited to 
different transversal states although in principle distinguishable interact statis- 
tically. It is an example of mutual statistical interaction [^2f introduced first 
to describe interlayer phase correlations in the double-layer Hall effect. The 
plateaux pattern for our anyonic vortices is given by (|3^). In distinction to 
bosonic Hall effect the constant a is nonzero. 

8 Summary 

We have given two new arguments why vortices in superconducting films should 
be anyons. In addition we have discussed two different experiments where this 
theoretical prediction could be verified. 
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FIGURE CAPTION. Detection of the Hall plateau. A$ is the moduli of 
the change in the external magnetic flux necessary to drive the change in the 
number of vortices by one. A$o is its average value far from the plateau, n is 
the density of vortices and no its value at the plateau. The curve " a" coresponds 
to decreasing while "b" to increasing number of vortices. The curve "c" is the 
difference between "a" and "b". 
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